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1 . INTRODUCTION 


The  work  reported  here  is  part  of  a continuing  program  aimed  at  deter- 
mining the  unsteady  flow  in  a transonic  compressor  stage.  The  stage  is 
installed  in  the  Turbopropulsion  Laboratories  of  the  Department  of  Aeronautics, 
Naval  Postgraduate  School. 

This  report  has  been  compiled  to  facilitate  use  of  the  data  acquisition 
and  analysis  programs  which  have  been  developed  primarily  for  the  study  of 
unsteady  fluctuating  pressures  on  the  casing  inner  wall.  The  equipment  and 
instrumentation  are  discussed  briefly  in  section  2.  The  data  acquisition 
system  and  programs  are  outlined  in  section  3 and  described  in  detail  in  the 
appendices.  The  post-real  time  data  analysis  programs  are  outlined  with  sample 
results  in  section  4 and  are  described  in  detail  in  the  appendices.  Section 
5 contains  conclusions  and  recommendations  for  further  work. 

2.  EQUIPMENT  AND  INSTRUMENTATION 
2.1  The  Transonic  Compressor 

The  transonic  compressor  test  rig  comprises  an  air  turbine  drive  unit 
and  an  induction  section  which  contains  a filter,  throttle,  settling  chamber 
and  flow  measuring  nozzle.  The  turbine  drive  unit  supplies  450  HP  at 
30,000  RPM.  The  compressor  is  designed  to  operate  at  30,460  RPM  with  a rela- 
tive tip  Mach  number  of  1.5.  At  the  design  RPM  and  the  tip  Mach  number,  the 
flow  angle  is  65®  and  the  pressure  ratio  is  1.6  at  a referred  flow  rate  of 
19  lbm/sec.  The  laboratory  facilities  and  the  test  rig  are  described  in 
detail  by  VAVRA  and  SHREEVE  (1972)  and  VAVRA  (1973). 
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2.2  Pressure  Measurement 

Eight  Kulite  CQL-080-25  pressure  transducers  with  natural  frequency  about 
125  kHz  are  mounted  with  their  diaphragms  flush  with  the  inner  case  wall  of 
the  compressor.  Further  details  are  reported  by  PAIGE  (1976).  Table  1 in 
Appendix  A gives  the  axial  and  circumferential  location  of  the  transducers 
relative  to  transducer  number  K6  which  is  the  furthest  upstream.  The  trans- 
ducers are  used  in  conjunction  with  Datel  Model  201C  instrumentation 
amplifiers  which  have  a flat  frequency  response  to  100  kHz. 

Each  Kulite  pressure  transducer  is  matched  by  a pneumatic  static  pressure 
tap  at  the  same  axial  location  in  the  case  wall  (except  in  one  case  - see  Table 
1 in  Appendix  A)  but  displaced  circumferentially.  Other  pneumatic  static  and 
total  pressure  taps  are  available  upstream.  A data  recording  system  (VAVRA 
and  SHREEVE,  1972)  is  used  to  record  both  the  steady  pressures  from  the 
pneumatic  taps  and  the  temperature  data.  The  paper  tape  output  from  this 
system  is  processed  using  a Hewlett  Packard  Model  HP9830A  programmable  calcu- 
lator to  provide  input  data  for  the  measurement  of  fluctuating  pressures  and 
to  establish  the  compressor  operating  point. 

2.3  The  Timing  Disk 

To  enable  synchronization  of  the  sampling  of  the  pressure  transducer 
outputs  with  the  rotation  of  the  rotor,  an  instrumented  timing  disk  is  fitted 
to  the  rotor  shaft.  The  disk  contains  holes  at  intervals  of  one  per  rotor 
blade  and  one  per  rotor  revolution.  Light  sensitive  diodes  and  wave  shaper 
circuits  provide  pulse  trains  to  control  sampling  of  the  pressure  transducers. 
This  system  is  described  in  detail  by  WEST  (1976). 


3.  DATA  ACQUISITION 
3.1  The  System  Hardware 

Figure  1 is  a schematic  of  the  data  acquisition  hardware  with  arrows 
indicating  the  flow  of  data  and  control  signals.  The  system  is  under  the 
control  of  the  HP  21MX  computer  which  operates  either  directly  or  through  the 
device  called  "Pacer"  to  control  the  analog-to-digital  (A/D)  converter 
(model  HP5610A)  and  which  transfers  data  to  the  HP9867B  mass  memory  unit  via 
the  HP9830A  calculator. 


The  peripheral  device  called  the  "pacer"  is  described  in  detail  by 
WEST  (1976)  who  originally  called  it  RPACE.  The  pacer  can  trigger  data 


acquisition  from  a stationary  transducer  at  any  fixed  point  in  the  rotating 
rotor  frame,  independent  of  the  rotor  speed.  In  effect,  it  divides  the 
circumference  of  the  rotor  into  9 intervals,  each  with  a circumferential 
length  equal  to  that  of  the  arc  (measured  at  the  wall)  across  two  adjacent 
blade  passages.  Each  of  the  9 intervals  are  subdivided  into  256  equal 
sub-intervals . 

The  pacer  receives  one  per  revolution  and  one  per  blade  input  signals 
from  the  timing  disk  and  performs  two  functions;  it  controls  the  timing  for 
data  acquisition  and  determines  the  speed  of  the  rotor. 

3.2  The  Program  KULITE 

KULITE  is  the  data  acquisition  program  (in  BASIC  language)  for  the 
HP  21MX  computer.  A flow  diagram,  listing,  variable  assignment  and  notes 
are  presented  in  Appendix  B.  KULITE  can  be  operated  in  three  modes,  viz. 
free-run,  calibration  and  pacer. 

In  free-run  mode  the  A/D  converter  operates  in  mode  4 (see  Reference 
No.  5).  Up  to  1616  samples  of  one  A/D  converter  channel  are  taken  with  a 
frequency  of  105  samples  per  second.  The  sampling  process  is  not  synchronized 
with  the  rotor  rotation. 

In  calibration  mode  the  A/D  converter  operates  in  mode  4 (see  Reference 
No.  5).  It  scans  through  A/D  converter  channel  numbers  1 to  12,  taking 
1616  samples  on  each  channel.  The  average  of  the  1616  samples  is  computed 
before  the  next  channel  is  sampled.  The  scan  is  performed  four  times,  with 
four  different  calibration  pressures  applied  to  the  reference  side  of  each 
Kulite  transducer. 
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In  pacer  mode  the  A/D  converter  operates  in  mode  0 (see  Reference  No. 


5).  Sampling  of  a Kulite  transducer  output  is  synchronized  with  rotor 
rotation  by  means  of  two  pulse  trains  generated  from  light  beams  chopped  by 
the  timing  disk  on  the  compressor  shaft.  One  pulse  train  has  a frequency  of 
one  per  rotor  revolution  and  the  other  has  a frequency  of  18  per  rotor 
revolution;  each  pulse  in  the  latter  train  corresponding  to  the  passing  of  a 
blade  past  a fixed  point.  A full  description  is  given  by  WEST  (1976). 

In  this  mode  a pressure  transducer  is  sampled  on  successive  revolutions 
at  a fixed  point  in  the  rotating  rotor  frame.  Currently,  the  sample  interval 
is  several  revolutions  of  the  rotor.  Changes  in  program  RPACE  would  allow 
samples  to  be  taken  at  intervals  of  one  revolution.  If  the  flow  can  be 
regarded  as  steady  in  the  rotating  rotor  frame  this  technique  enables 
measurement  of  the  wall  pressure  distribution  "carried  around"  by  the  rotor. 
Flow  unsteadiness  in  the  rotating  frame  can  be  averaged  or  the  frequency  con- 
tent of  the  unsteadiness  in  successive  samples  can  be  examined.  In  this 
report  only  averaged  data  from  10  samples  taken  at  each  of  128  points  across 
two  rotor  blade  passages,  is  presented. 

In  all  three  modes  of  operation  the  program  KULITE  transfers  data  from 
the  HP21MX  to  the  HP9830A. 

3.3  The  Program  TRAN4 

TRAN4  is  the  data  acquisition  program  (in  BASIC  language)  for  the  HP9830A 
programmable  calculator.  It  receives  data  from  the  HP21MX  computer,  processes 
it  and  stores  data  on  a disk  of  the  HP9867B  mass  memory.  A flow  diagram, 
listing,  variable  assignment  and  notes  are  presented  in  Appendix  C. 
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3.4  The  Program  RF.SET1 

RESET1  initializes  a record  number  on  the  storage  disk  so  that  at  the 
start  of  a run  data  can  be  stored  in  file  DATAY1  beginning  at  the  first  record. 
The  program  is  listed  in  Appendix  D. 

4.  DATA  ANALYSIS 
4.1  The  Program 

Off  line  data  analysis  is  at  present  performed  on  the  HP9830A  with  the 
BASIC  language  programs  MAPI,  MAP2 , CONT,  CONTI,  PLOTSA,  PLOTSB  and  TITIPK. 
These  programs  are  described  in  detail,  with  listings,  flow  diagrams  and 
notes,  in  Appendices  E through  J. 

MAPI  is  used  to  determine  the  sensitivity  of  the  Kulite  pressure  trans- 
ducers from  data  acquired  with  KULITE  in  the  calibration  mode.  In  addition, 
MAPI  is  used  to  convert  the  voltages  sampled  at  the  pressure  transducer  out- 
puts to  pressure  coefficients. 

MAP2  is  used  to  convert  the  8 x 128  array  of  measured  pressure  coeffi- 
cients to  a 29  x 128  array  through  quadratic  interpolation  in  the  axial 
direction.  The  program  was  written  to  reduce  the  effects  of  the  course  trans- 
ducer spacing  in  the  axial  direction.  However,  care  must  be  exercised  when 
it  is  used  to  interpolated  across  discontinuities  such  as  shock  waves  and 
rotor  blades.  Linear  interpolation  is  available  through  use  of  the  program 
CONT  or  CONTI  to  plot  contours  of  casing  wall  pressures. 

CONT  Is  used  to  plot  contours  of  constant  casing  wall  pressure  (in  the 
frame  of  the  rotor)  from  an  array  of  pressure  coefficients.  (i.e.  it  pro- 
duces a wall  pressure  "map") . The  program  will  accept  any  general  rectangular 
array  provided  that  the  spacing  in  each  direction  is  uniform.  This  latter 
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requirement  restricts  it's  use  in  this  application  to  arrays  obtained  from  MAP2. 

CONTI  is  used  to  plot  contours  of  constant  casing  wall  pressure  when  the 
array  of  measured  pressure  coefficients  contains  nonuniform  spacing  in  the 
axial  direction.  Nonuniform  spacing  results  in  this  application  from  the 
axial  location  of  the  Kulite  pressure  transducers. 

PLOTSA  is  used  to  plot  (on  the  HP9862A  plotter)  the  uncalibrated  pressure 
distribution  (in  volts)  across  a blade  pair  for  a given  Kulite  transducer. 

The  input  data  is  that  originating  from  pacer  mode  of  operation.  The  program 
is  also  used  to  plot  the  output  of  the  one  per  blade  signal  from  the  timing 
disk. 

PLOTSB  is  used  to  plot  (on  the  HP9862A  plotter)  the  uncalibrated  free- 
run  data  (in  volts)  from  a given  transducer  against  circumferential  distance. 

TITIPK  is  used  to  superimpose  the  blade  tip  profiles  on  the  wall 
pressure  maps. 

4.2  Sample  Results 

The  results  presented  here  are  intended  only  to  illustrate  the  capabili- 
ties of  the  programs.  Comprehensive  results  will  be  given  in  a subsequent 
report. 

Figure  2 is  a plot  versus  circumferential  distance  of  the  average 
pressure  in  the  frame  of  the  rotor  across  an  arc  of  the  casing  wall  equivalent 
to  two  blade  passages.  The  pressure  coefficient  is  defined  in  Appendix  E. 

The  plot  was  made  with  PLOTSA  using  data  acquired  in  the  pacer  mode  of  operation. 
The  precise  location  of  each  distribution  relative  to  the  rotor  blades  is  not 
defined  here.  The  locations  are  known  approximately  from  the  blade  pair 
number  specified  in  the  acquisition  program.  They  are  located  precisely  from 
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Circumferential  distance  in  opposite  direction  to  rotation 


Figure  2 

I 


. Waveshapes  of  unsteady  pressure  distributions  across  tvo  blade 
passages  (uncalibrated  and  with  arbitrary  offsets)  for  the  Kulite 
transducers.  Data  taken  in  pacer  mode.  502  design  speed;  throttled 
to  near  surge.  8.7  lbm/sec  referred  flow  rate.  Pressure  ratio  - 
1.155:1.  Blade  pair  # 2. 
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Figure  3.  Typical  record  of  raw  one  per  blade  signal  taken  in  pacer 
mode.  Blade  passage  frequency  is  4.55  kHz. 


Figure  4.  Typical  record  of  data  taken  in  free-run  mode  from  Kulite  number  K10. 
The  record  is  comprised  of  1616  samples  taken  at  a frequency  of 
100  kHz.  Each  cycle  is  due  to  a blade  passage  with  a blade  passing 
frequency  of  4.55  kHz.  Compressor  operating  conditions  as  in  figure  2 


Circumferential  (opposite  to  direction  of  rotation) 


K6  K7  K8  K9  KIO  Kl!  K 1 2 KI3 


^ Streamwise  Axial 

Figure  5.  Contours  of  constant  pressure  coefficient  ACp  plotted 

by  CONTI.  (See  Appendix  E for  definition  of  AC  ) . Blade  pair  number 
2.  Compressor  operating  conditions  as  in  figuri  2. 
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Figure  6.  Smoothed  contours  of  constant  pressure  coefficient  obtained  from  Figure  5. 
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Figure  7.  Contours  of  constant  pressure  coefficient  plotted  by  CONT  with 
25x255  array.  Blade  pair  number  2.  Compressor  operating 
conditions  as  in  Figure  2. 
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Figure  8.  Contours  of  constant  pressure  coefficient  plotted  by 

CONT  with  7x64  array  as  subset  of  array  used  for  Figure  7. 
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Circumferential  (opposite  to 


Figure  9.  Contours  of  constant  pressure  coefficient  plotted  by  CONT  with 

8x128  array.  Note  there  is  a small  error  in  blade  location.  Blade 
pair  number  2.  Sixty  percent  design  speed.  Throttled  to  near  surge. 
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a knowledge  of  the  orientation  of  the  timing  disk  relative  to  the  rotor  in 
conjunction  with  the  phase  relationship  between  the  one  per  blade  timing 
signal  and  the  pressure  distributions  from  the  transducers.  A typical  record 
of  the  one  per  blade  signal,  taken  directly  from  the  photo-diode  output  in 
the  pacer  mode  is  plotted  in  figure  3. 

The  distributions  shown  are  of  output  voltage  from  the  Kulite  transducers 
which  at  this  stage  have  not  been  scaled  to  take  account  of  the  tranducer 
calibrations.  The  rapid  changes  in  the  signals,  in  particular  those  from 
transducers  K9,  K10  and  Kll,  are  due  to  the  passage  of  a rotor  blade  across  a 
transducer  and  provide  one  means  of  estimating  blade  location.  The  distribu- 
tions have  been  plotted  by  linear  interpolation  between  the  128  points  across 
the  blade  pair.  Each  of  the  128  points  has  been  obtained  from  the  average  of 
ten  samples,  one  sample  being  taken  approximately  each  tenth  revolution  of 
the  rotor. 

Figure  4 is  a typical  record  of  data  taken  in  the  free-run  mode  from  a 
particular  transducer.  The  fundamental  frequency  is  that  of  the  blade 
passage  past  the  transducer  (4.55kHz).  The  plot  is  a linear  interpolation 
between  1616  samples  taken  at  a frequency  of  100  kHz. 

Figure  5 is  a representative  map  of  the  contours  of  constant  pressure 
coefficient  with  respect  to  an  upstream  reference  static  pneumatic  pressure. 

It  was  generated  with  the  program  CONTI  from  the  8 x 128  array  of  measurements 
which  were  taken  in  the  pacer  mode  and  calibrated  with  the  program  MAPI.  The 
map  across  two  adjacent  rotor  blade  passages  is  thus  that  of  the  mean  wall 
pressures  "carried  around"  by  the  rotor.  The  method  of  calculating  the 
the  pressure  coefficients  from  the  Kulite  transducer  outputs  and  the  measure- 
ments from  the  pneumatic  taps  is  given  in  Appendix  E. 
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The  location  of  the  blade  tip  in  figure  5 was  determined  from  the  circum- 
ferential location  of  the  rapid  change  in  pressure  distribution  associated 
with  transducer  K9  (figure  2).  The  transducer  K9  has  an  axial  location 
such  that  it  is  crossed  by  the  leading  edge  of  rotor  blades.  A more  precise 
location  of  blades  can  be  obtained  from  measurement  of  the  circumferential 
location  of  the  one  per  blade  raw  signal  (figure  3)  relative  to  the  pressure 
distributions  (figure  2)  and  from  knowledge  of  the  location  of  the  one  per 
blade  holes  in  the  timing  disk  relative  to  the  rotor. 

Smoothing  of  the  Wall  Pressure  Maps 

It  is  clear  in  figure  5 that  linear  interpolation  between  pressure 
coefficients  in  the  axial  direction  across  a blade  gives  incorrect  contours 
along  the  blade  between  the  transducer  locations.  This  difficulty  can  only 
be  overcome  satisfactorily  by  using  more  Kulite  transducers  to  provide  a 
finer  mesh  than  the  rather  coarse  one  provided  by  eight  transducers.  In  this 
study  the  contours  in  the  vicinity  of  the  blades  have  been  smoothed  graphically 
by  connecting  with  smooth  curves  those  points  in  the  map  at  which  the  pressure 
coefficient  is  the  same  and  which  lie  on  the  lines  scanned  by  the  transducers. 
At  operating  conditions  which  give  rise  to  shock  waves  it  is  possible  that  a 
similar  smoothing  procedure  will  be  required.  In  the  long  run  there  is  a need 
for  a numerical  interpolation  technique  which  avoids  interpolation  across 
blades. 

Figure  7 is  a wall  pressure  map  produced  by  CONT  for  similar  operating 
conditions  to  those  used  to  obtain  figure  5.  However,  pressures  were 
measured  at  255  points  in  the  circumferential  direction  and  MAP2  was  used  to 
generate  pressure  coefficients  at  25  equally  spaced  stations  in  the  axial 
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direction  between  transducers  K6  and  K12.  It  is  clear  that  quadratic  inter- 
polation does  not  avoid  the  problems  of  contour  distortion  by  interpolation 
across  a blade.  In  view  of  the  very  slow  execution  of  CONT  with  a 25  x 255 
array  of  pressure  coefficients,  it  is  recommended  that  an  8 x 128  array  be 
used  with  CONTI. 

Figure  8 is  a wall  pressure  map  of  the  same  data  as  in  figure  7 except 

that  only  every  fourth  circumferential  point  (of  the  original  255)  was  used 

in  the  reduction  with  CONTI  (i.e.  axial  interpolation  is  linear  through  use 

of  CONTI  but  MAP2  is  not  used).  The  resulting  7 x 64  array  of  pressure 

coefficients  gives  rise  to  very  similar  contours  to  those  shown  in  figure  7 

for  the  finer  mesh.  In  fact,  the  appearance  of  contours  with  C = 0.4  near 

P 

the  blades  in  figure  7 might  be  an  erroneous  result  of  quadratic  interpolation 
across  blades.  Note  that  the  grid  lines  in  figure  8 which  indicate  transducer 
locations  are  incorrectly  plotted  to  have  equal  spacing,  but  this  does  not 
affect  the  above  comments. 

Figure  9 is  a plot  of  contours  for  a higher  compressor  rotational  speed. 
The  data  reduction  techniques  were  identical  to  those  used  to  obtain  figure  5. 
Note  that  the  blade  location  shown  was  that  calculated  for  the  data  of  figure 
5.  The  location  is  slightly  in  error  due  to  a difference  in  the  tuning  of  the 
phase-lock  loop  circuit  in  the  pacer  for  the  second  test.  This  problem  has 
been  solved  recently  by  using  the  near-discontinuity  in  the  pressure  distribu- 
tion measured  by  the  transducer  at  the  leading  edge  to  position  the  blade  tip. 

Again  it  is  stressed  that  the  results  presented  are  prelinimary  and  are 
intended  merely  to  demonstrate  the  methods  used  and  the  capabilities  of  the 
system. 


20 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Programs  for  the  acquisition  and  reduction  of  fluctuating  casing  wall 
pressures  in  a transonic  compressor  stage  have  been  developed  and  run 
successfully  with  the  compressor  at  this  time  operating  at  up  to  60  percent 
of  the  design  speed.  Evaluation  of  the  data  acquisition  system  on  a 
mechanical  simulator  indicates  that  it  can  operate  over  the  full  speed 
range  of  the  compressor.  In  fact,  at  higher  speeds  the  signal  to  noise 
ratio  in  the  Kulite  transducer  outputs  will  improve  significantly  because  of 
the  higher  pressures  that  will  be  encountered. 

The  pacer  mode  of  synchronized  sampling  has  made  it  possible  to  deter- 
mine in  a versatile  manner  the  wall  pressure  maps  in  the  rotating  frame  of 
the  rotor.  The  pacer  system  can  also  be  used  to  obtain  measurements  of 
flow  properties  away  from  the  wall,  e.g.  flow  velocity  measurements  with  a 
dynamic  probe.  Wall  pressure  maps  have  been  presented  solely  to  demonstrate 
the  capabilities  of  the  pacer  technique  and  of  the  data  acquisition  and 
analysis.  Comprehensive  data  will  be  presented  and  interpreted  in  a 
subsequent  report. 

There  are  some  aspects  of  the  programs  which  can  be  refined  or  which 
need  further  evaluation  and  the  following  recommendations  are  made. 

1.  The  subroutine  RPACE  in  KULITE  causes  a sample  to  be  taken  in  pacer  mode 
about  every  tenth  revolution  of  the  rotor  at  60  percent  of  design  speed. 
This  causes  a delay  in  data  acquisition  which  could  be  reduced  by 
modifying  subroutine  RPACE. 

2.  The  degree  of  steadiness  of  pressure  distributions  in  the  frame  of  the 
rotor  needs  further  investigation.  This  should  begin  with  an  examination 
of  the  standard  deviations  already  computed  in  pacer  mode.  In  separate 
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tests,  a larger  number  of  samples  (at  least  500)  should  be  taken  for  each 


i 

of  several  steps  between  blade  pairs  and  the  variations  at  each  step 
examined  for  frequency  content.  The  measurements  should  then  be  repeated 
with  the  case  wall  rotated  peripherally  by  at  least  90°. 

3.  The  technique  of  calibration  of  the  Kulite  pressure  transducers  under 
operating  conditions  effectively  takes  account  of  change  in  transducer 
sensitivity  with  temperature.  Change  in  transducer  offset  (d.c.  level) 
with  temperature  is  presently  handled  by  equating  the  time-average  trans- 
ducer output  voltage  with  the  steady  pressure  obtained  from  a pneumatic 
static  tap  at  the  same  axial  location.  The  relationship  between  the 
steady  pressure  indicated  by  the  pneumatic  tap  and  the  time  time-averaged 
pressure  at  the  tap  needs  further  investigation. 

4.  Because  transducer  K10  and  pneumatic  tap  S10  are  not  coincident  axially 
it  is  necessary  to  interpolate  between  readings  at  S10  and  Sll.  The 
interpolation  in  MAPI  is  at  present  linear  but  its  adequacy  has  not  been 
fully  evaluated. 

5.  The  large  pressure  gradients  in  the  axial  direction  across  rotor  blades 
are  not  resolved  well  because  of  the  limited  number  of  transducers.  Two 
additional  transducers,  located  midway  between  K9  and  K10  and  K10  and 
Kll,  would  greatly  alleviate  this  difficulty.  Linear  axial  interpolation 
across  blades,  as  in  CONTI,  is  misleading  and  hand  smoothing  of  contours 
near  blades  is  presently  necessary.  Quadratic  interpolation,  as  in  MAP2 , 
does  r.ot  solve  the  problem.  Extrapolation  of  data  up  to  but  not  across 

a blade  surface  should  be  investigated. 

j 
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6.  Shock  waves  have  not  been  encountered  at  the  low  operating  speeds  at 
which  the  present  data  was  obtained.  The  accuracy  of  resolution  of 
shock  waves  should  be  studied  in  the  light  of  the  above  discussion  re- 
garding large  pressure  gradients  across  blades. 

7.  The  blades  have  been  located  on  the  wall  pressure  maps  in  this  report 
from  knowledge  of  the  point  in  the  circumferential  pressure  distribution 
(indicated  by  transducer  K9  at  the  blade  leading  edge)  at  which  the 
circumferential  pressure  gradient  is  steepest.  This  technique  is  subject 
to  an,  as  yet,  undetermined  uncertainty  due  to  irregularities  in  the 
geometry  from  blade  to  blade.  The  alternative  procedure,  whereby 
blades  are  located  by  use  of  the  phase  relationship  between  the  one  per 
blade  signal  and  the  circumferential  pressure  distributions,  also  needs 
further  evaluation. 

8.  In  its  present  form  the  data  acquisition  system  requires  frequent  key- 
board entries  by  the  operator.  In  principle  the  system  can  be  fully 
automated  by  pre-entering  all  necessary  data  with  DATA  statements  and  by 
replacing  INPUT  statements  by  READ  statements.  Some  WAIT  statements  in 
KULITE  would  be  needed  to  allow  TRAN4  to  catch  up  to  KULITE. 

9.  The  format  of  graphical  outputs  can  be  improved  by  using  the  plotter 
to  add  alphameric  information. 

10.  Two-way  data  transfer  beteen  the  HP21MX  and  the  HP9830A  is  feasible. 

This  capability  should  be  developed  to  enable  use  of  the  faster  HP21MX 
for  repetitious  data  reduction. 
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APPENDIX  A 

Table  1 contains  the  axial  and  circumferential  location  of  the  Kulite 
pressure  transducers  and  the  axial  location  of  the  pneumatic  static  pressure 


Kulite 

transducer 

number 


Pneumatic 
static  tap 
numb  er 


Axial  distance 
downstream  of  K6 
(inches) 


Circumferential 
location  relative 
to  K6  in 

direction  of  rotation 


- 

0 

- 

0.50 

- 

1.00 

- 

1.37 

— 

1.75 

— 

2.12 

- 

2.50 

- 

3.00 

S6 

0 

S7 

0.50 

S8 

1.00 

S9 

1.37 

S10 

1.55 

Sll 

2.12 

S12 

2.50 

S13 

3.00 

Table  1 Location  of  Kulite  pressure  transducers  and  pneumatic  static  pressure 
taps. 
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APPENDIX  B 
DETAILS  OF  KULITE 

KULITE  Is  the  data  acquisition  program  for  the  HP21MX  computer.  Its  three 
modes  of  operation  are  indicated  in  section  3.2.  Figure  10  is  a flow  diagram 
of  the  program  and  a listing  is  given  in  Table  2. 


Variable  Assignment  for  KULITE 


11 

12 

13 

14 

15 

16 
17 
13 
19 


A1 

T1 

N1 

N2 

M 

S 

R 

L 

A3, A4 ,A5, A6 
A[  101,  16] 
B[101,  16] 

C [ 10] , D[ 10 ] 


- Run  It.  Same  as  Run  It  in  Log  Book  assigned  to  each 
start-up  of  the  compressor 

- Test  It.  Refers  to  a particular  operating  condition 
within  a run. 

- Day 

- Month 

- Year 

-A/D  converter  mode  It. 

- Samples /channel  in  free- run  mode. 

- Not  used. 

- Experiment  It.  Refers  to  either 

(i)  One  time  series  of  free-run  data, 

(ii)  A complete  set  of  calibration  readings  (averaged) 
for  all  transducers,  or 

(iii)  Averaged  pressures  across  one  blade  pair  in  Pacer 
mode . 

- Channel  It.  Refers  to  A/D  converter. 

- Transducer  It. 

- Samples /point  in  Pacer  mode. 

- Blade  pair  It. 

- Mean  of  pressure  samples  at  a point  in  Pacer  mode. 

- Standard  deviation  of  pressure  samples  at  a point  in 
Pacer  mode . 

- Mean  of  one/blade  signal  at  a point  in  Pacer  mode. 

- Row  number  in  K matrix  of  calibrations. 

- Associated  with  subroutine  RPACE  and  defined  by  WEST  (1976) 

- Consecutive  free-run  samples  stored  row  by  row. 

- Buffer  in  subroutine  R5610 

- Buffers  in  subroutine  R5610 
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E[4,  255] 


K[5,  12] 


This  is  the  matrix  of  averaged  data  taken  in  RPACE 
across  two  blade  passages.  Note  that  128  points  are 
taken  across  two  blade  passages.  This  number  can  be 
changed  to  255  by  changing  Line  195. 

- Matrix  of  calibrations.  Rows  1,  2,  3,  each  contain 
averaged  calibration  voltages  for  the  twelve  trans- 
ducers. Each  of  rows  1,  2,  3 corresponds  to  a 
different  calibration  level.  Row  5 contains  II,  12, 

13,  14,  15,  19  and  the  three  reference  pressures (which 
are  keyed  in  on  request)  in  K [ 5 , 7 ] , K[5,8],  K [ 5 , 9 ] . Row 
4 is  treated  as  another  calibration  level  and  is  used 
to  scan  the  offsets  if  needed.  In  that  case  any  value 
can  be  input  to  K[5,10]  for  P Ref. 
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Table  2. 


Listing  of  KULITE 


I nEM  "KULITE"  DATA  LOGGING  PROGRAM-- SIMMONS--  11  APRIL,  1977 
5 DI v At  I 0 1,  1 60,  or  I 0 1,  1 61  , Ct  I 0]  , DC  101,  EC  4,  2550, XC  5/  1 20 

9 PRINT  "ENTER RUN#, MONTH, DAY, YEAR" 

10  FOR  1=1  TP  5 

I I FPP  J= 1 TO  12 

12  LET  KC  I,  J0  = 0 

13  NEXT  J 

14  NEXT  I 

15  INPUT  11,14,13,15 

16  PD I NT  "ENTER  0 IF  9830  IS  RESET 1 FUR  END" 

17  INPUT  R8 

18  IF  R8>.5  THEN  47  0 

20  PRINT  "ENTER TEST#,  EXPERIMENT#" 

25  INPUT  12,19 

26  PnINT  "0  FOR  CALIBRATION OTHER’*  I SE  1" 

27  INPUT  Z 1 

28  PRINT#  8J  Z 1 

20  I F Z l<  . 5 THEN  1 29 

30  PPINT  "ENTER MURE#,  CHANNEL#,  TRANSDUCER#" 

35  INPUT  I 6,  A 1 , T 1 
40  PPINT#  8J  I I 
45  PRINT#  8.'  I 2 
50  PRINT#  8;  I 3 
55  PRINT#  8J  14 
60  PRINT#  8,'  IS 
65  PRINT#  81  I 6 
70  PPINT#  8,'  19 
75  PP INT#  8 J T 1 
PC  °PINT#  5*A| 

85  IF  16=0  THEN  165 

90  IF  16=4  THEN  103 

95  PRINT  "VPPMG  MODE#" 

98  GOTO  3 0 

103  PnINT  "ENTER  SAMPLES  / CHANNEL  (NOT  > 1616)" 

10  5 INPUT  17 
110  PCINT#  8J  17 

115  R5  61 0(7, At  1,  11,  17, Al,  16, BC  1,  11) 

116  °P  I NT  "DATA  TAKEN" 

120  FO®  J* 1 Tn  101 

124  FOP  I»1  TO  16 

125  PRINT#  8,*ACJ,n 

126  NEXT  I 

127  NEXT  J 

128  GOTO  16 

129  LFT  L=  0 

130  LET  L*L*  1 

131  P°INT  "ENTER  P REF (IN.  H20  REL  ATMOS)  IF  °EADY  FOP  CALIBRATION- 

132  I NPVT  Kt  5, 6+L 1 
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Table  2.  Cont 


133  FOR  J 1=  1 TP  12 

134  PRINT  " START  IMG  CALIBRATION  OF  A/D  CHANNEL  * "J  I 
13  5 L FT  A 1 3 J 1 

13  6 P5  6 1 0 (7,  AC  1 , 11,  I 6 1 6,  Al,  4,  3t  1,  11) 

137  LET  XtL,  J 11*0 

138  FPP  J* 1 TP  101 

139  TCP  1=1  TO  15  STEP  2 

140  IF  ABSfAt J,  I 1 >< .98  THEN  145 

141  PRINT  -********************•*******" 

142  PC!NT  “*  * ** ** OVERLOAD  ON  "J  1 "******■ 

143  PRINT  -************************•***“ 

144  GOTO  131 

145  LET  XtL, J ll*XtL, J Il+AtJ, I) 

146  NEXT  I 

147  NEXT  J 

148  LET  XtL,  J ll=XtL,  Jll/1616 

149  NEXT  J 1 

150  P°INT  "0  FOP  ANOTHER  CAL  I BRATI  ON- -OTHER’*’  I SE  1* 

151  INPUT  Z 2 

152  IF  Z2<  . 5 THEN  1 30 

153  LET  Xt  5,  11*1  1 

154  LET  XT  5.21=12 

155  LET  Xt  5,  31  = 13 

156  LET  Xt  5,  41  = I 4 

157  LET  Xt  5. 51* I 5 
153  LET  Xt  5.  61*  19 

159  FOP.  1=1  TO  5 

160  FOR  J* 1 TO  12 

161  popiTi  8 j Kt  I , J 1 

162  NEXT  J 

163  NEXT  I 

164  GOTO  16 

165  PRINT  "ENTER  BLADE  PAIR  *,  SANPLES/PC INT" 

166  INPUT  N2,Nl 

167  PRINT*  8JN1 

168  PRINT*  3JN2 
170  LET  A3*  0 

180  IF  Tl»INT(Tl/2)*2*. 1 THEN  190 

131  LET  A6*32768*N 2*256-64 

182  GOTO  195 

190  LET  A6*32763*N 2*256 
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Table  2.  Cont. 


195  F0C  1= ! TO  255  STEP  2 
20  0 LET  A3  = A6+  I 
215  LET  P=0 
225  RPACEC  A3^A4/A5) 

230  R5610C7,  CC  M,NI,  Al,  0,  DC  M> 

240  FOR  J=1  TO  Ml 
2M  LET  9C  J,  1 1 = CCJ] 

242  NEXT  J 

243  LET  A3=32768*N2*256+I 

244  RPACEC  A3.»A4*A5) 

245  R5610(7,CC  1 LML  Oi  0<  DC  11) 

246  FOR  J * 1 TO  Ml 

247  LET  RaR*CCJ) 

250  NEXT  J 

25  5 LET  RaR/Nl 
260  LET  S=0 
270  LET  V=0 
2R  0 FOP  J = 1 TO  M 1 
290  LET  N = M*BC  J,  1 1 
300  NEXT  J 
310  LET  M*M/N  1 
320  FOR  J=  1 TO  N 1 

330  LET  S*S-K CECJ, 11-M)*CBCJ, 11-M) ) 
340  NEXT  .J 

350  LET  S* SQR( S/ ( N 1 - 1 ) > 

360  LET  EC  1,  Il-M 
370  LET  EC  2,  I1*S 
380  LET  EC  3,  I1  = A4 

38  5 LET  EC  4,  I )*R 

39  0 NEXT  I 

40  0 FOR  J*  1 TO  4 

41 0 FOR  1=1  TO  255  STEP  2 

420  PRINT#  8; EC J, II 

43  0 NEXT  I 

440  NEXT  J 

450  GOTO  16 

47  0 END 
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Notes  on  KULITE 


1.  A/D  converter  channels.  It  is  essential  that  the  "raw"  one  per  blade 
signal  be  input  to  channel  1 of  the  A/D  converter.  Allocation  of  the 
other  channels  is  not  unique  but  the  allocation  in  Table  3 is  recommended. 
Channels  11  and  12  are  scanned  but  at  present  are  not  used  in  subsequent 
analysis. 


A/D  Converter 
Channel  Number 

Signal 

0 

one  per  blade  raw  signal 

1 

K6 

2 

K7 

3 

K8 

4 

K9 

5 

K10 

6 

Kll 

7 

K12 

8 

K13 

* 

9 

P c - P 

ret  atmos 

10 

S2  - P r 
ref 

11 

Unused 

12 

Unused 

Table  3.  Allocation  of  signals  to  A/D  converter  channels. 


P ^ is  pressure  applied  to  reference  side  of  Kulite  transducers. 

2.  Subroutine  R5610  is  described  by  WEST  (1976,  p.  17). 

3.  In  calibration  mode  the  scan  through  the  twelve  channels  must  be  made  four 
times.  The  first  scan  must  be  with. the  pressure  tapping  S2  applied 
simultaneously  to  the  reference  ride  of  Kulite  transducers.  The  second 


32 


and  third  scans  must  be  made  with  other  steady  pressures  applied  to  the 
reference  side  of  the  Kulites.  The  fourth  scan  must  be  made  to  satisfy 
the  program  but  at  this  stage  the  data  taken  is  not  used  in  subsequent 
analysis.  This  scan  is  included  to  enable  logging  of  the  offsets  on  the 
Kulite  amplifiers  should  they  be  of  interest. 

The  program  searches  for  overloads  (i.e.  greater  than  0.98  volts  or  less 
than  - 0.98  volts  in  the  calibration  signals^  If  it  detects  an  overload 
among  alternate  samples  in  the  1616  samples  taken  from  any  transducer  the 
offending  A/D  converter  channel  number  is  displayed,  the  scan  is  aborted 
and  the  program  is  reset  to  repeat  the  scan.  The  limit  of  0.98  volts  can 
be  changed  in  line  140. 

The  one  per  revolution  signal  from  the  timing  wheel  indicates  the  origin 
for  circumferential  measurements  around  the  rotor.  The  pacer  then  uses 
the  one  per  blade  signal  to  divide  the  rotor  circumference  into  9 equal 
intervals,  the  first  interval  beginning  at  the  origin.  These  intervals 
are  designated  by  blade  pair  numbers,  although  the  start  of  an  interval 
need  not  coincide  with  a blade  tip.  Each  interval  represents  a circum- 
ferential length.  In  the  rotor  frame,  equal  to  that  of  the  arc  (measured 
at  the  wall)  across  two  adjacent  blade  passages.  Each  of  the  nine 
intervals  is  divided  into  256  sub- intervals . In  pacer  mode  the  scan 
across  an  interval  begins  after  the  first  sub-interval  and  ends  after  the 
255th  sub-interval.  With  stepping  sequentially  across  the  sub- intervals 
in  pairs,  a total  of  128  points  are  sampled.  It  is  convenient  to  take  10 
samples  at  each  point  (one  sample  approximately  each  ten  revolutions)  to 
compute  the  mean  and  standard  deviation. 


6.  Even  numbered  transducers  (e.g.  K6,  K8  etc.)  are  located  on  one  axial 
line  and  odd  numbered  transducers  are  on  another  axial  line  which  is  dis- 
placed around  the  casing  wall  by  10  degrees  in  the  direction  opposite 
that  of  rotation  of  the  rotor.  The  parity  of  the  transducer  number  is 
evaluated  in  line  180. 

The  variable  A3  determines  the  time  (in  terms  of  degrees  of  rotation 
of  the  rotor)  after  the  one  per  revolution  pulse  when  a sample  is  taken 
at  point  I.  For  example,  in  line  190,  A6  = 32768  + N2  * 256  for  odd 
transducer  numbers,  and  A3  = A6  + I . This  defines  the  sampling  time 
(approximately  each  tenth  revolution)  for  point  I in  blade  pair  N2. 

(I  = 1 to  255,  N2  = 1 to  9)  . Point  I can  be  sampled  10  degrees  earlier 
for  even  numbered  transducers  by  setting 
A6  = 32768  + N2  * 256  - 64 
and  A3  = A6  + I 

The  subroutine  RPACE  is  described  in  more  detail  by  WEST  (1976)  . 

7.  During  each  scan  of  a transducer  across  a blade  pair  the  raw  one  per  blade 
is  also  sampled.  This  is  used  later  in  TRAN4  to  determine  the  location 

of  the  measured  pressure  distribution  relative  to  the  rotor. 
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APPENDIX  C 


DETAILS  OF  TRAN4 


TRAN4  is  the  data  acquisition  program  for  the  HP  9830A  programmable  calculator. 


Figure  11  is  a flow  diagram  and  a listing  is  in  Table  4. 


Variable  Assignment  for  TRAN4 


N1 

Z1 

A3 


RECORD  //.  i.e.  Number  of  first  available  record  in  DATAY1 
IDENTIFIER  ( * 0 FOR  CALIBRATION  - OTHERWISE  1) 

BLADE  PAIR  It 


A2 


SAMPLES /POINT  in  Pacer  mode. 


T1 

T2 

A[5, 128] 


LOCATION  (between  1 and  255)  of  point  in  blade  pair  where 
1/blade  signal  is  0.5  volts  and  increasing. 


Location  of  point  in  blade  pair  where  1/blade  signal  is  0.5 
volts  and  decreasing. 


M. 


step  2 


255 


255 


A4 


255 


1 2 255 

II,  12,  13,  14,  15,  16,  17,  Tl,  Al,  A2,  A3,  0,  0,  0 


NOTE  that  variable  names  in  this  array  are  those  used  in  21MX 
program.  They  should  not  be  confused  with  variables  used  in 
TRAN4.  This  is  the  matrix  of  averaged  data  taken  in  RPACE 
across  two  blade  passages.  Note  that  128  points  can  be  changed 
by  changin  dimension  statement  and  the  FOR  loop. 

B [ 9 ] - Buffer  for  identification  data 

[II,  12,  13,  14,  15,  16,  19,  Tl,  A1]T 

NOTE  that  the  variable  names  in  this  array  are  those  used  in 
21MX  program. 
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D [ 16 , 102] 


K[5,  12] 

RECY  # 
DATAY1 


Free  run  data 
Stored  row  by  row 


I 


I 


11 

12 

13 

14 

15 

16 
19 
T1 
A1 
17 


Column  102 


NOTE  that  variable  names  in  column  102  are  those  used  in 
21KX  program. 

Same  as  K[5,  12]  in  21MX  program. 


DATA  FILES 


This  file  contains  1 record.  It  contains  N1  which  is  the 
number  of  the  1st  available  record  in  file  DATAY1 

Data  file  for  calibration,  free  run  and  Pacer  data.  Le.  for 
K,  D and  A arrays  respectively.  It  contains  300  records. 
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Table  4.  Listing  of  TKAN4 
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1 .0  t op  i i rr  • 
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i--j  next  i 

r O OEM  7[  ,T  MODI  4 
’■  r,u  [F  Bit-']- it  rnt(1  s -;ij 
■V  lr  £:[■  j - 4 THEN  3*50 
30  DISP  "MROl  j (10 BE  # " ; 

; ,.i  -EH  - c.'  E I E *-sr-»r»  T 0 F • . C r'tl.  1 8K  i 
.00  I- OR  I - ! TO  =, 

. 10  t'tjP  .i  ! TO  1. 
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. 40  1 j t j T 0 t.y 
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■ •:  0 - N TtF  'll*  I o 
; ’0  t-'OR  I = i T 1 - 

. 0 K 0 F J - i T 0 1 O 1 

: enter  i . - rr  i - j ] 

4M0  NEXT  i 
410  NEXT  I 
4:0  F 0 R 1 = 1 T O 0 
4 '0  ii[  I . 104  ]--£:[  I ] 

440  NE  :T  I 
4 50  £' [ 10)  1 0 P 1=17 
460  •-EmB  tell*  Ml 
4~0  HRT  PRINT  te  J:?D 
400  N 1 =N  1 -t- 1 3 

4 00  REmB  # 1 » 1 

500  PRINT  41  I HI 
f,  0 , 0 T 0 050 

500  REM  RECEIVE  RftCEP  DhTP 

5 SO  ENTER  1 1 « ♦ MX 

540  ENTER  ■ 1 < ■ m 

5t:0  FOR  J=1  TO  4 
5X0  FOR  1=1  TQ  K.3 
5“0  ENTER  • 1 - ’ .'Ml  J>  I J 
5 aO  NEXT  I 

5 R0  NEXT  j 

t.00  REIT  FIND  POOH  ION  OF  RErEPEMOE 
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Table  4.  Cont. 
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BEST  AVAILABLE  COPY 


Notes  on  TRAN4 


1. 


2. 


(line  30)  DATAY1  is  a file  on  a removable  disk  for  temporary  storage.  At 
the  end  of  Run  number  n (n  is  a two  digit  integer)  a file  CKRWn  must 
be  opened  on  the  fixed  disk  and  data  must  be  copied  into  it  from  DAYAl 
for  long  term  storage.  The  number  of  records  in  CKRWn  must  equal  the  sura 
of  record  number  printed  out  with  identification  of  the  last  experiment 
and  a number  k where 

k = 1 if  last  experiment  was  a calibration 

=5  if  last  experiment  was  in  pacer  mode 

= 13  if  last  experiment  was  in  free-run  mode 
Figure  12  contains  sample  print-out  from  TRAN4  for  the  three  modes. 

(Line  600)  The  position  of  the  centre  of  the  raw  one  per  blade  pulse  is 
found  relative  to  the  pressure  distribution  across  a blade  pair  by 
searching  through  the  128  averaged  samples  at  each  point  to  find  the 
sample  which  first  exceeds  0.5  volts  and  the  first  subsequent  sample 
which  is  less  than  0.5  volts.  The  corresponding  point  numbers  are 
averaged.  If  the  pulse  is  not  found  (due  to  inadequate  signal  level) , 
the  program  displays  an  ERROR.  After  the  correct  one  per  blade  signal 
has  been  re-established  both  KULITE  and  TRAN4  must  be  rerun. 
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APPENDIX  L) 


DETAILS  OF  RESET1 

RESET1  initializes  to  1 the  number  in  the  file  named  RECY//.  This 
enables  data  acquired  at  the  start  of  a run  to  be  stored  at  the  start  of  the 
file  named  DATAY1 . Table  5 is  a listing. 


Table  5.  Listing  of  RESET1 


cE.i:  i • : I -,'IHOM 
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* • r N f ri 


AVA  i 

i i k i 
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APPENDIX  E 


DETAILS  OF  MAPI 

MAPI  is  used  to  compute  the  sensitivities  of  the  Kulite  pressure  trans- 
ducers and  to  convert  voltages  sampled  at  the  pressure  transducer  outputs 
to  pressure  coefficients.  Figure  13  is  a flow  diagram  and  a listing  is  in 
Table  6. 


Variable  Assignment  for  MAPI 


N 

T 

PI 

P2 

M 

U1 

A1 

Q 

Cl,  C2,  C3 

RO 

TO 

A 

PO 


compressor  RPM 
total  temperature  (called  T. 


TOT 


elsewhere)  measured  at  axial 


location  of  S2  (entered  in  degrees  F) . 

Static  pressure  measured  at  S2  (entered  in  inches  of 

water  absolute) . 

Total  pressure  (pT0T)  measured  at  axial  location  of  S2 
(entered  in  inches  of  water  absolute) . 

square  of  Mach  number  at  axial  location  of  S2. 

2 2 

square  of  rotor  tip  speed  in  ft  /sec  . 

square  of  speed  of  sound  at  axial  location  of  S2  in  ft/sec. 

Reference  pressure  for  computing  pressure  coefficients. 

First,  second  and  third  calibration  pressures  applied  to 
reference  side  of  Kulite  transducers  (inches  of  water 
relative  to  atmospheric) . 

Record  number  for  calibration  data  on  disk. 

Kulite  transducer  number. 

Kulite  transducer  sensitivity  in  inches  of  water  per  volt. 

Pressure  from  pneumatic  wall  static  tap  corresponding  in  axial 
location  to  a Kulite  ansducer.  (inches  of  water  relative 
to  pressure  at  S2) . 


Figure  13.  Flow  diagram  for  MAP  1 

44 


Table  6.  Listing  of  MAPI 
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Table  6.  Cont 


R1 


N 0 


U 

L 

AS[5,  128]  - 

C[5,  12] 

P [ 32  ] 

U [8 , 4] 

L[8,  4] 

DATAY1 


Record  number  for  pacer  data. 

N0+  1 is  pacer  data  quarter  number.  The  128  samples  are 
divided  into  4 sets  of  32. 

Maximum  pressure  coefficient  in  set  across  two  blade  passages. 
Minimum  pressure  coefficient  in  set  across  two  blade  passages. 
Same  as  A[5,  128]  in  TRAN4 . 

Same  as  K[5,  12]  in  TRAN4 . 

Array  of  pressure  coefficients  for  one  transducer  and  in  one 
quarter  of  pressure  distribution  across  two  blade  passages. 

U[I,  J]  is  local  maximum  pressure  coefficient  for  transducer 
number  I and  quarter  J. 

L[I,  J]  is  local  minimum  pressure  coefficient  for  transducer 
number  I and  quarter  J. 

Same  data  file  as  in  TRAN4. 


TRESS  - contains  pressure  coefficients.  The  32  records  contain  in 

order  quarters  1,  2,  3,  4 for  transducer  6,  quarters  1,  2, 
3,  4 for  transducer  7,  etc. 


Notes  on  MAPI 

1.  The  correct  CKRWn  file  name  for  the  run  under  consideration  must  be 
entered. 

2.  Storage  of  pressure  coefficients.  The  128  pressure  coefficients  asso- 
ciated with  each  transducer  have  been  grouped  into  4 sets  (quarters)  with 
each  set  being  stored  in  a separate  record.  This  has  been  done  to 
facilitate  use  of  interpolation  programs  such  as  MAP2.  Interpolation 
expands  the  size  of  the  data  array  so  that  pressure  coefficients  must 

be  recalled  from  the  mass  memory  in  subsets  in  order  to  meet  the  storage 
limitations  of  the  HP9830A. 
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3.  Sensitivities  of  the  Kullte  transducers  are  computed  as  follows: 


Cl,  C2,  C3  are  the  three  steady  pressures  applied  to  the  reference  side 
of  the  Kulite  transducers.  They  correspond  to  mean  output  voltages 
C[l,  TO  - 5],  C[2,  TO  - 5],  C [ 3 , TO  - 5]  from  Kulite  transducer  number 
TO.  The  sensitivity  of  a transducer  (i.e.  the  slope  of  its  calibration 
curve  at  a particular  mean  operating  temperature  and  pressure)  is  ob- 
tained from  a least  squares  fit  of  a straight  line  through  the  three 
points . 

4 . Calculation  of  pressure  coefficients . 

The  Kulite  transducer  output  voltages  E are  converted  to  pressure 

coefficients  C as  follows: 

P 

C = ((E  - E)  * A + P0)/Q 
P 

where  E = transducer  mean  output  voltage  obtained  during  calibration 
with  S2  on  reference  side  of  diaphragm. 

A = sensitivity  of  transducer  in  inches  of  water  per  volt. 

PO  =*  mean  wall  pressure  (at  same  axial  location  as  Kulite  trans- 
ducer) measured  with  pneumatic  tap.  (inches  of  water  relative 
to  S2) . 

Q = reference  dynamic  pressure  (inches  of  water  absolute)  computed 
as  follows: 

The  reference  dynamic  pressure  is  expressed  in  terms  of  the  upstream 
density  p and  the  upstream  flow  velocity  measured  in  the  rotating  frame 
of  the  rotor.  Hence 


1 2 2 

q*P  (v4  + in 

where  V » flow  velocity  at  station  S2  (ft/sec). 
U * rotor  tip  speed 


N 


nN  5.5 
30  x 12 

rotor  RPM 


ft. sec. 


(1) 
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It  follows  (noting  that  variable  names  are  not  necessarily  the  same  as 
in  the  listing  of  MAPI)  that 


1 2 2 2 
Q - ~ y?  (M  + U /a  ) 


(2) 


where  P,  a,  M are  static  pressure,  speed  of  sound  and  Mach  number  at 
station  S2. 


But 


(3) 


and 


(4) 


where  Xj,  and  are  total  temperature  and  pressure  at  S2. 

By  introducing  (1),  (3)  and  (4)  into  (2),  Q can  be  calculated  in  terms 

of  P,  P_  and  T_  . 

, T T 

5.  Note  that  S10  and  K10  are  not  at  the  same  axial  location.  For  purposes 
of  computing  pressure  coefficients  an  effective  mean  wall  pressure  at  K.10 
is  obtained  by  linear  interpolation  using  values  at  S10  and  Sll . 

6.  The  maximum  and  minimum  pressure  coefficients  are  computed  to  aid  in 
the  choice  of  contours  when  using  the  programs  CONT  or  CONTI. 


t 
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APPENDIX  F 


DETAILS  OF  MAP2 


MAP2  is  used  to  convert  the  8 x 128  array  of  measured  pressure  coeffi- 
cients to  a 29  x 128  array  through  quadratic  interpolation  in  the  axial 
direction.  A listing  is  in  Table  7. 

Variable  Assignment  for  MAP2 


B [ 32  ] 

temporary  storage  of  pressure  coefficients 

C (29 , 32  ] 

array  of  interpolated  pressure  coefficients  across 
quarter  of  a blade  pair. 

one 

P [8 , 32] 

Array  of  pressure  coefficients  at  Kulite  transducer 
across  one  quarter  of  blade  pair. 

locations 

X [ 8 ] 

Axial  location  of  Kulite  transducers  downstream  of 
ducer  K6.  (inches). 

trans- 

Q 

quarter  number. 

PRESS 

Same  data  file  as  in  MAPI. 

INTER 

File  for  storage  of  interpolated  pressure  coefficients  across 
one  quarter  of  blade  pair  (15  records). 

Notes  on  MAP2 

1.  Lagrangian  interpolation  is  used,  i.e.  if  P^,  P0,  and  P^  are  known  at 


x^,  x^,  x^  then 


( x—  x 2 ) (x-x^) 


( x— x, ) (x-x_) 

P(x)  - 7 ^-7-  P + *—r~. P 

(x. -x„) (x, -x,)  1 (x_-x, ) (x--x,)  2 


(x-x.)  (x-x.) 

— p 

(x3~x1) (x3~x2)  3 


Table  7-  Listing  of  MAP 2 


i 
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* 


.‘•i  run 


2.  As  a result  of  the  equispaced  interpolation,  the  transducer  axial  loca- 
tions will  in  general  not  coincide  with  axial  stations  in  the  interpolated 

array. 

3.  Three  interpolations  are  made  in  the  axial  direction  between  transducer 
measurements . 

4.  The  program  must  be  run  for  each  quarter  of  the  array  of  measured 
pressure  coefficients.  After  each  running  of  the  program  the  contours 
must  be  plotted  with  CONT  prior  to  running  MAP2  for  another  quarter. 

5.  Interpolation  between  X[l]  and  X [ 2 ] is  made  with  a quadratic  through 
pressure  coefficients  at  X[l],  X [ 2 ] and  X[3J.  Interpolation  between 
X [ 7 ] and  X [ 8 ] is  made  with  a quadratic  through  X[6],  X [ 7 ] and  X [ 8 ] . 
Interpolation  between  X[I]  and  X [ I + 1]  (for  I > 1 and  < 6)  is  made  by 
averaging  the  quadratic  through  X[I  - 1],  X[I]  and  X[I  + 1]  and  the 
quadratic  through  X[I],  X[I  + 1]  and  X[  1+2]. 
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APPENDIX  G 
DETAILS  OF  CONTI 

CONTI  is  used  to  plot  contours  of  constant  casing  wall  pressure  from 
an  array  of  measured  pressure  coefficients.  The  program  handles  the  non- 
uniform  axial  spacing  of  the  Kulite  transducers  but  requires  uniform 
circumferential  spacing  in  the  array.  CONTI  is  written  to  accept  the  array 
generated  by  MAPI.  Figure  14  is  a flow  diagram  and  a listing  is  in  Table  8. 

Variable  Assignment  for  CONTI 

P[I,  J]  - array  of  pressure  coefficients 

A[3,  3],  B[3,  3],  Z[3],  Q [ 3]  , D{2,  2],  F[2],  G[2,  2],  H[2]  - defined  in 

note  6 of  this  appendix. 

axial  location  (inches)  of  Kulite  transducers  downstream 
of  K6 . 

dimensions  of  P[I,  J]  in  axial  and  circumferential  directions 
respectively. 

value  of  pressure  coefficient  on  a contour, 
triangular  element  number. 

number  of  starting  element  in  a contour  plot, 
number  of  finishing  element  in  a contour  plot. 

Same  file  as  in  MAPI. 


1.  The  X (axial)  and  Y (circumferential)  dimensions  of  P must  be  entered 
before  program  is  run.  P[8,  128]  is  nearly  the  maximum  array  size  that 
can  be  stored  in  the  HP9830A. 

2.  The  axes  are  drawn  so  that  X (axial)  runs  from  K6  to  K13  and  Y (circum- 
ferential) runs  from  the  start  to  the  end  of  a blade  pair.  Contours 
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Notes  on  CONTI 
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will  not  quite  go  to  these  extremities  because  the  pacer  system  starts 
sampling  at  1/256  th  of  a blade  pair  and  finishes  sampling  at  255/256th 
of  a blade  pair.  Recall  that  the  start  of  a blade  pair  need  not  coin- 
cide with  a blade  because  of  phase  lags  in  the  pacer  system  and  the 
location  of  the  timing  disk  relative  to  the  rotor.  Grid  lines  parallel 
to  the  Y-axis  are  the  lines  scanned  by  the  Kulite  transducers. 

3.  If  other  transducer  locations  are  used  their  axial  distance  downstream 
of  K6  must  be  entered  before  the  program  is  run. 

4.  The  variable  plotted  in  the  Y-direction  is  the  column  number  in  P.  The 
variable  plotted  in  the  X-direction  is  the  axial  location  of  a transducer 
(downstream  of  K6)  and  is  derived  from  the  corresponding  row  number  in  P. 
When  setting  up  the  X-Y  plotter  it  is  advisable  to  make  both  the  X and  Y 
scales  equal  to  twice  full  scale.  Note  that  the  circumferential  distance 
along  the  wall  across  two  blade  passages  is  3.8A7  inches. 

5.  Entering  of  the  first  and  last  element  numbers  enables  faster  plotting 
of  a contour  which  is  known  in  advance  to  cover  only  a limited  part  of 
the  field. 

6.  Triangular  element  representation  of  the  surface  defining  the  pressure 
coefficient  distribution,  Cp  (X,Y) 

The  surface  - Cp  (X,Y)  is  approximated  by  triangular  elements  as 

illustrated  in  figure  15.  Element  numbers  E are  as  shown.  This  process 

represents  linear  interpolation  between  measured  pressure  coefficients. 

Contours  are  obtained  from  the  intersection  of  planes  C =*  constant  with 

P 

this  approximation  to  the  pressure  distribution.  Thus  the  contours  are 
composed  of  straight  line  segments.  The  nodes  correspond  in  the  Y-directicn 
(circumferential)  to  points  at  which  pacer  data  is  available.  The  node 
numbering,  in  the  local  sense,  for  typical  odd  and  even  numbered  elements 
is  shown  in  figure  15. 
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Stceam^ija  ax<al  X 


Figure  L5  Numbering  of  triangular  elements  in  CONT  and  CONTI.  Also 

numbering  of  nodes  (vertices  of  triangles)  for  add  and  even 
numbered  elements. 
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If  Z = aX  + bY  + c is  the  plane  containing  the  three  nodes  of  a 


triangle  then 


Z(l) 

XI 

Yl 

1 

Z ( 2 ) 

= 

X2 

Y2 

1 

Z(31 

X3 

Y3 

X. 

where  Z(I)  is  the  pressure  coefficient  at  local  node  number  I (i,e.  at 
XI,  YI). 

This  matrix  equation  is  written  as 

Z = Aa 

and  hence 

0 - A ^Z  = BZ  say. 


If  E is  odd,  the  node  coordinates  are  related  to  the  element  number 
as  follows: 

XI  - A[l,  1]  = INT  ( ( (E  + l)/2  - .00001)/(J0  - 1) ) + 1 
Yl  - (E  + 1 ) / 2 - (A[l,  1]  -1)  * (JO  -1) 

X2  - XI  , Y2  - Yl  + 1 
X3  - XI  + 1 , Y3  » Yl 

Note  that  the  number  .00001  is  included  to  avoid  problems  associated  with 
round-off  error. 

Z - aX  + bY  + c is  plane  containing  triangle.  For  Z = c',  line  of 
intersection  with  triangle  (i.e.  contour)  is  given  by  aX  + bY  + (c  -c')  = 
On  the  side  X = I » A[l,  I],  (See  Figure  15) 

Y - (-  aA[l,  1]  - (c  - c'))/b 

- (-  Q [ 1 ] * A[l,  1]  - Q ( 3 ) + c')/Q(2) 
provided  Q(2)  i 0 
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If  Q(2)  = 0 , the  contour  is  parallel  to  the  side  in  question.  On  the 

side  Y - J - A [ 1 , 2] 

X - (-Q [2 ] * A [1 , 2]  - Q ( 3)  + c')/Q(l) 
provided  Q ( 1 ) i4  0 

On  the  hypotenuse 

X + Y = A [ 3 , 1]  + A[3,  2] 
aX  + bY  - c'  - c 
Hence 


1 

i 

X 

m 

A[ 3,  1]  + A[ 3,  2] 

a 

b. 

.Y. 

c1  - c 

This  matrix  equation  is  written  as 
DH  = F 

Hence  H = GF  where  G = D ^ 

am 

If  E is  even,  the  node  coordinates  are  related  to  the  element  number 
as  follows: 

XI  = A[l,  1]  - INI  ((E/2  - .00001) /(JO  - 1))  + 1 

Y1  =■  A [ 1 , 2]  - E/2  - (A[l,  1]  - 1)  * (JO  - 1)  + 1 

X2  - XI  + 1 , Y2  - Y1  - 1 

X3  - XI  + 1 , Y3  - Y1 

Z = aX  + bY  + c is  plane  containing  the  triangle.  For  Z * c'  , line  of 

intersection  (i.e.  contour)  is  given  by  aX  + bY  + (c  - c')  = 0 . 

On  the  side  X * I * A[2,  1]  , 

Y - (-  Q(l)  * A [2 , 1]  - Q ( 3)  + c’)/Q(2) 
provided  Q ( 2)  + 0 

On  the  side  Y - J - A[l,  2] 

X - (-  Q(2)  * A [1,  2]  - Q( 3)  + c')/Q(l) 

provided  Q(l)  t 0 . 
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. 


On  the  hypotenuse 

X+  Y » A[2,  1]  + A[2,  2] 
aX  + bY  = c ' - c 


Hence 


1 1 

a b 


LYJ 


A[2,  1]  + A [2 , 2] 


c - c 


This  matrix  equation  is  written  as 


DH  = F 

Hence  H = GF  where  G = D 

The  intersections  of  the  contour  with  the  sides  of  the  triangular 
elements  are  computed  in  the  manner  outlined  above.  Tests  are  performed 
to  identify  those  points  of  intersection  which  lie  on  side  of  the  triangle 
as  opposed  to  intersections  which  lie  on  extrapolations  of  the  sides 
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APPENDIX  H 


DETAILS  OF  CONT 

CONT  is  used  to  plot  contours  of  constant  casing  wall  pressure  from 
an  array  of  pressure  coefficients  in  which  both  the  axial  and  the  circum- 
ferential spacings  are  uniform.  It  can  be  used  to  plot  contours  for  an 
array  output  by  MAP2.  The  program  differs  from  CONTI  only  in  its  plotting 
of  x-coordinates  of  contours.  A listing  is  in  Table  9. 

Notes  on  CONT 

1.  If  MAP2  is  used  for  interpolation,  the  file  PRESS  must  be  replaced  bv 
file  INTER  and  P must  have  dimensions  P[29,  32]  . Contours  must  be 
plotted  by  processing  the  array  of  interpolated  pressure  coefficients 
in  four  sections.  This  requires  shifting  of  origin  on  the  x-y  plotter 
for  each  quarter. 

2.  The  variables  plotted  in  the  x and  y directions  are  respectively 

the  row  number  and  column  number  of  the  sub-array  of  pressure  coef ficients. 
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Table  9 


Listing  of  CONT 
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APPENDIX  I 

DETAILS  OF  PLOTS  A AND  PLOTS  3 

PLOTS  A and  PLOTS  B are  used  to  plot  pacer  raw  data  and  free-run 
raw  data  respectively  against  circumferential  distance.  The  programs  are 
listed  in  Tables  10  and  11.  Note  that  PLOTS  A,  with  I * 1 in  line  110, 
plots  the  output  of  the  designated  Kuli'a  transducer.  If  I = 4 the  raw 
one  per  blade  signal  on  n/ D converter  channel  number  0 is  plotted. 
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Table  11.  Listing  of  PLOTSB 
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APPENDIX  J 


DETAILS  OF  TITIPK 

TITIPK  is  used  to  superimpose  hlade  tip  profiles  on  wall  pressure 
contours.  The  program  is  loaded  onto  the  programmable  keys  of  the  HP9830A. 
The  program  is  listed  in  Table  12.  The  blade  tip  profile  is  tabulated  by 
PAIGE  (1976),  figure  2 and  table  1).  The  axial  and  circumferential  units 
are  inches  and  the  program  is  compatible  with  a wall  pressure  map  which  has 
dimensions  of  3 inches  axially  by  3.847  inches  circumferentially.  The 
"lower  left"  and  "upper  right"  on  the  plotter  should  be  set  to  the  corre- 
sponding points  on  the  wall  pressure  map. 

The  key  programs  should  be  "continued"  after  CONTI  has  been  run  so 
that  the  Kulite  data  is  available  in  main  memory.  When  <CONT>  <f^>  is 
issued,  the  location  (YO)  of  the  blade  leading  edge  from  the  lower  boundary 
of  the  wall  maps  is  calculated  from  the  data  of  transducer  K9  and  appears  in 
the  display.  <CONT>  <f^>  *s  then  issued  and  YO  is  requested  as  an  input. 

The  blade  profiles  are  then  drawn  and  the  key  program  ends.  The  contour 
plotting  with  CONTI  can  be  continued  by  issuing  <CONT>  1250. 
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APPENDIX  K 


OPERATING  PROCEDURE  FOR  DATA  ACQUISITION  SYSTEM 

1.  Load  Real  Time  Executive  Basic  into  HP21MX  computer. 

2.  Load  KULITE  into  HP21MX.  Tune  pacer. 

3.  Put  disk  labeled  "Transonic  Compressor  - Paige"  into  HP9867B  mass  memory. 

4.  At  start  of  a run,  get  RESET1  from  unit  0 and  run  RESET1  to  initialize 
the  number  in  file  R£CY It  to  unity. 

5.  Scratch  RESET1. 

6.  Get  TRAN4  from  unit  0. 

7.  Run  TRAN4.  The  HP9830A  display  will  remain  blank  while  the  HP9830A 
awaits  data  from  the  HP21MX. 

8.  Run  KULITE,  noting  that  the  two  mode  switches  on  the  Pacer  must  be  set 
according  to  the  A/D  converter  mode  to  be  used  (i.e.  0 or  4) . In  cali- 
bration, four  scans  of  the  twelve  channels  must  be  made.  The  first  scan 
must  be  with  the  pressure  at  S2  on  the  reference  side  of  the  Kulite 
transducers.  The  second  and  third  scans  are  made  with  other  steady 
calibration  pressures  on  the  reference  side.  The  fourth  scans  must  be 
made  but  any  signals  can  be  used  on  the  A/D  converter  provided  that  they 
do  not  cause  overloads.  In  mode  0 the  pacer  must  not  be  altered  during 
an  experiment. 

9.  On  completion  of  an  experiment  both  the  HP9830A  and  the  HP21MX  must  be 
reset  as  instructed  by  their  displays  prior  to  performing  another 
experiment. 

10.  On  completion  of  a run  the  data  which  is  stored  temporarily  in  DATAY1 
file  must  be  transferred  to  permanent  files  CKRWm  duplicated  on  both 
unit  0 and  unit  1.  First  open  CRKWm  where  m is  the  run  number.  The 


73 


length  of  CKRWm  must  be  set  at  k records  where 
k » record  number  for  last  experiment  + g 


and  g * 1 if  last  experiment  is  a calibration 

5 if  last  experiment  is  a pacer  experiment 

13  if  last  experiment  is  a free-run  experiment 

11.  To  abort  the  HP21MX  program,  enter  AB. 
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